Truncating variants of the sterol recognition region of SHH cause hypertelorism phenotype rather than hypotelorism-holoprosencephaly.

Sonic hedgehog signaling molecule (SHH) is a key molecule in the cilia-mediated signaling pathway and a critical morphogen in embryogenesis. The association between loss-of-function variants of SHH and holoprosencephaly is well established. In mice experiments, reduced or increased signaling of SHH have been shown to be associated with narrowing or excessive expansion of the facial midline, respectively. Herein, we report two unrelated patients with de novo truncating variants of SHH presenting with hypertelorism rather than hypotelorism. The first patient was a 13-year-old girl. Her facial features included hypertelorism, strabismus, telecanthus, malocclusion, frontal bossing, and wide widow's peak. She had borderline developmental delay and agenesis of the corpus callosum. She had a nonsense variant of SHH: Chr7(GRCh38):g.155802987C > T, NM_000193.4:c.1302G > A, p.(Trp434*). The second patient was a 25-year-old girl. Her facial features included hypertelorism and wide widow's peak. She had developmental delay and agenesis of the corpus callosum. She had a frameshift variant of SHH: Chr7(GRCh38):g.155803072_155803074delCGGinsT, NM_000193.4:c.1215_1217delCCGinsA, p.(Asp405Glufs*92). The hypertelorism phenotype contrasts sharply with the prototypical hypotelorism-holoprosencephaly phenotype associated with loss-of-function of SHH. We concluded that a subset of truncating variants of SHH could be associated with hypertelorism rather than hypotelorism.

Pathogenicity evaluation of variants of uncertain significance at exon-intron junction by splicing assay in patients with Mowat-Wilson syndrome.

High-throughput sequencing has identified vast numbers of variants in genetic disorders. However, the significance of variants at the exon-intron junction remains controversial. Even though most cases of Mowat-Wilson syndrome (MOWS) are caused by heterozygous loss-of-function variants in ZEB2, the pathogenicity of variants at exon-intron junction is often indeterminable. We identified four intronic variants in 5/173 patients with clinical suspicion for MOWS, and evaluated their pathogenicity by in vitro analyses. The minigene analysis showed that c.73+2T>G caused most of the transcripts skipping exon 2, while c.916+6T>G led to partial skipping of exon 7. No splicing abnormalities were detected in both c.917-21T>C and c.3067+6A>T. The minigene analysis reproduced the splicing observed in the blood cells of the patient with c.73+2T>G. The degree of the exon skipping was concordant with the severity of MOWS; while the patient with c.73+2T>G was typical MOWS, the patient with c.916+6T>G showed milder phenotype which has been seldom reported. Our results demonstrate that mRNA splicing assays using the minigenes are valuable for determining the clinical significance of intronic variants in patients with not only MOWS but also other genetic diseases with splicing aberrations and may explain atypical or milder cases, such as the current patient.

Breakpoints in complex chromosomal rearrangements correspond to transposase-accessible regions of DNA from mature sperm.

Constitutional complex chromosomal rearrangements (CCRs) are rare cytogenetic aberrations arising in the germline via an unknown mechanism. Here we analyzed the breakpoint junctions of microscopically three-way or more complex translocations using comprehensive genomic and epigenomic analyses. All of these translocation junctions showed submicroscopic genomic complexity reminiscent of chromothripsis. The breakpoints were clustered within small genomic domains with junctions showing microhomology or microinsertions. Notably, all of the de novo cases were of paternal origin. The breakpoint distributions corresponded specifically to the ATAC-seq (assay for transposase-accessible chromatin with sequencing) read data peak of mature sperm and not to other chromatin markers or tissues. We propose that DNA breaks in CCRs may develop in an accessible region of densely packaged chromatin during post-meiotic spermiogenesis.

Impaired gating of γ- and ε-AChR respectively causes Escobar syndrome and fast-channel myasthenia.

OBJECTIVE: To dissect the kinetic defects of acetylcholine receptor (AChR) γ subunit variant in an incomplete form of the Escobar syndrome without pterygium and compare it with those of a variant of corresponding residue in the AChR ε subunit in a congenital myasthenic syndrome (CMS). METHODS: Whole exome sequencing, α-bungarotoxin binding assay, single channel patch-clamp recordings, and maximum likelihood analysis of channel kinetics. RESULTS: We identified compound heterozygous variants in AChR γ and ε subunits in three Escobar syndrome (1-3) and three CMS patients (4-6), respectively. Each Escobar syndrome patient carries γP121R along with γV221Afs*44 in patients 1 and 2, and γY63* in patient 3. Three CMS patients share εP121T along with εR20W, εG-8R, and εY15H in patients 4, 5, and 6, respectively. Surface expressions of γP121R- and εP121T-AChR were 80% and 138% of the corresponding wild-type AChR, whereas εR20W, εG-8R, and εY15H reduced receptor expression to 27%, 35%, and 30% of wild-type εAChR, respectively. γV221Afs*44 and γY63* are null variants. Thus, γP121R and εP121T determine the phenotype. γP121R and εP121T shorten channel opening burst duration to 28% and 18% of corresponding wild-type AChR by reducing the channel gating equilibrium constant 44- and 63-fold, respectively. INTERPRETATION: Similar impairment of channel gating efficiency of a corresponding P121 residue in the acetylcholine-binding site of the AChR γ and ε subunits causes Escobar syndrome without pterygium and fast-channel CMS, respectively, suggesting that therapy for the fast-channel CMS will benefit Escobar syndrome.

Clinical diversity and molecular mechanism of VPS35L-associated Ritscher-Schinzel syndrome.

PURPOSE: The Retriever subunit VPS35L is the third responsible gene for Ritscher-Schinzel syndrome (RSS) after WASHC5 and CCDC22. To date, only one pair of siblings have been reported and their condition was significantly more severe than typical RSS. This study aimed to understand the clinical spectrum and underlying molecular mechanism in VPS35L-associated RSS. METHODS: We report three new patients with biallelic VPS35L variants. Biochemical and cellular analyses were performed to elucidate disease aetiology. RESULTS: In addition to typical features of RSS, we confirmed hypercholesterolaemia, hypogammaglobulinaemia and intestinal lymphangiectasia as novel complications of VPS35L-associated RSS. The latter two complications as well as proteinuria have not been reported in patients with CCDC22 and WASHC5 variants. One patient showed a severe phenotype and the other two were milder. Cells established from patients with the milder phenotypes showed relatively higher VPS35L protein expression. Cellular analysis found VPS35L ablation decreased the cell surface level of lipoprotein receptor-related protein 1 and low-density lipoprotein receptor, resulting in reduced low-density lipoprotein cellular uptake. CONCLUSION: VPS35L-associated RSS is a distinct clinical entity with diverse phenotype and severity, with a possible molecular mechanism of hypercholesterolaemia. These findings provide new insight into the essential and distinctive role of Retriever in human development.

Angelman syndrome with mosaic paternal uniparental disomy suggestive of mitotic nondisjunction.

Angelman syndrome (AS) is caused by the functional absence of the maternal ubiquitin-protein ligase E3A (UBE3A) gene. Approximately 5% of AS is caused by paternal uniparental disomy of chromosome 15 (UPD(15)pat), most of which is considered to result from monosomy rescue. However, little attention has focused on how UPD(15)pat occurs. We suggest the mitotic nondisjunction mechanism as a cause of UPD(15)pat in a six-year-old patient presenting with distinctive characteristics in line with AS. DNA methylation screening of 15q11-q13 showed a paternal band and a faint maternal band, suggestive of mosaic status. By trio-based microsatellite analysis, we confirmed a large proportion of UPD(15)pat cells and a small proportion of cells of biparental origin. Single nucleotide polymorphism (SNP) microarray revealed isodisomy of the entire chromosome 15. These results suggest that the UPD(15)pat of the patient resulted from mitotic nondisjunction, which may also be the cause of other cases of AS with UPD(15)pat.

Atypical Sotos syndrome caused by a novel splice site variant.

Sotos syndrome is usually caused by haploinsufficiency of NSD1; it is characterized by overgrowth, craniofacial features, and learning disabilities. We describe a boy with Sotos syndrome caused by a splicing variant (c.4378+5G>A). The clinical manifestations included severe connective tissue involvement, including joint hypermobility, progressive scoliosis, pectus deformity, and skin hyperextensibility; no overgrowth was observed.

Whole-exome analysis of 177 pediatric patients with undiagnosed diseases.

Recently, whole-exome sequencing (WES) has been used for genetic diagnoses of patients who remain otherwise undiagnosed. WES was performed in 177 Japanese patients with undiagnosed conditions who were referred to the Tokai regional branch of the Initiative on Rare and Undiagnosed Diseases (IRUD) (TOKAI-IRUD). This study included only patients who had not previously received genome-wide testing. Review meetings with specialists in various medical fields were held to evaluate the genetic diagnosis in each case, which was based on the guidelines of the American College of Medical Genetics and Genomics. WES identified diagnostic single-nucleotide variants in 66 patients and copy number variants (CNVs) in 11 patients. Additionally, a patient was diagnosed with Angelman syndrome with a complex clinical phenotype upon detection of a paternally derived uniparental disomy (UPD) [upd(15)pat] wherein the patient carried a homozygous DUOX2 p.E520D variant in the UPD region. Functional analysis confirmed that this DUOX2 variant was a loss-of-function missense substitution and the primary cause of congenital hypothyroidism. A significantly higher proportion of genetic diagnoses was achieved compared to previous reports (44%, 78/177 vs. 24-35%, respectively), probably due to detailed discussions and the higher rate of CNV detection.

Target enrichment long-read sequencing with adaptive sampling can determine the structure of the small supernumerary marker chromosomes.

Structural analysis of small supernumerary marker chromosomes (sSMCs) has revealed that many have complex structures. Structural analysis of sSMCs by whole genome sequencing using short-read sequencers is challenging however because most present with a low level of mosaicism and consist of a small region of the involved chromosome. In this present study, we applied adaptive sampling using nanopore long-read sequencing technology to enrich the target region and thereby attempted to determine the structure of two sSMCs with complex structural rearrangements previously revealed by cytogenetic microarray. In adaptive sampling, simple specification of the target region in the FASTA file enables to identify whether or not the sequencing DNA is included in the target, thus promoting efficient long-read sequencing. To evaluate the target enrichment efficiency, we performed conventional pair-end short-read sequencing in parallel. Sequencing with adaptive sampling achieved a target enrichment at about a 11.0- to 11.5-fold higher coverage rate than conventional pair-end sequencing. This enabled us to quickly identify all breakpoint junctions and determine the exact sSMC structure as a ring chromosome. In addition to the microhomology and microinsertion at the junctions, we identified inverted repeat structure in both sSMCs, suggesting the common generation mechanism involving replication impairment. Adaptive sampling is thus an easy and beneficial method of determining the structures of complex chromosomal rearrangements.

The novel and recurrent variants in exon 31 of CREBBP in Japanese patients with Menke-Hennekam syndrome.

Menke-Hennekam syndrome-1 (MKHK1) is a congenital disorder caused by the heterozygous variants in exon 30 or 31 of CREBBP (CREB binding protein) gene mapped on 16p13.3. It is characterized by psychomotor delay, variable impairment of intellectual disability (ID), feeding difficulty, autistic behavior, hearing impairment, short stature, microcephaly, and facial dysmorphisms. The CREBBP loss-of-function variants cause Rubinstein-Taybi syndrome-1 (RSTS1). The function of CREBBP leading to MKHK1 has not been clarified so far, and the phenotype of MKHK1 significantly differs from that of RSTS1. We examined six patients with de novo pathogenic variants affecting the last exon of CREBBP, and they shared the clinical features of MKHK1. This study revealed that one frameshift and three nonsense variants of CREBBP cause MKHK1, and inferred that the nonsense variants of the last exon could further help in the elucidation of the etiology of MKHK1.

Integrated diagnosis based on transcriptome analysis in suspected pediatric sarcomas.

Pediatric solid tumors are a heterogeneous group of neoplasms with over 100 subtypes. Clinical and histopathological diagnosis remains challenging due to the overlapping morphological and immunohistochemical findings and the presence of atypical cases. To evaluate the potential utility of including RNA-sequencing (RNA-seq) in the diagnostic process, we performed RNA-seq in 47 patients with suspected pediatric sarcomas. Histopathologists specialized in pediatric cancer re-evaluated pathological specimens to reach a consensus diagnosis; 42 patients were diagnosed with known subtypes of solid tumors whereas 5 patients were diagnosed with undifferentiated sarcoma. RNA-seq analysis confirmed and refined consensus diagnoses and further identified diagnostic genetic variants in four of the five patients with undifferentiated sarcoma. Genetic lesions were detected in 23 patients, including the novel SMARCA4-THOP1 fusion gene and 22 conventional or recently reported genetic events. Unsupervised clustering analysis of the RNA-seq data identified a distinct cluster defined by the overexpression of rhabdomyosarcoma-associated genes including MYOG and CHRNG. These findings suggest that RNA-seq-based genetic analysis may aid in the diagnosis of suspected pediatric sarcomas, which would be useful for the development of stratified treatment strategies.

Establishment of in-hospital clinical network for patients with neurofibromatosis type 1 in Nagoya University Hospital.

Neurofibromatosis type 1 (NF1) is a genetic multisystem disorder. Clinicians must be aware of the diverse clinical features of this disorder in order to provide optimal care for it. We have set up an NF1 in-hospital medical care network of specialists regardless of patient age, launching a multidisciplinary approach to the disease for the first time in Japan. From January 2014 to December 2020, 246 patients were enrolled in the NF1 patient list and medical records. Mean age was 26.0 years ranging from 3 months to 80 years. The number of patients was higher as age at first visit was lower. There were 107 males (41%) and 139 females. After 2011, the number of patients has increased since the year when the medical care network was started. Regarding orthopedic signs, scoliosis was present in 60 cases (26%), and bone abnormalities in the upper arm, forearm, and tibia in 8 cases (3.5%). Neurofibromas other than cutaneous neurofibromas were present in 90 cases (39%), and MPNST in 17 cases (7.4%). We launched a multidisciplinary NF1 clinic system for the first time in Japan. For patients with NF1, which is a hereditary and systemic disease associated with a high incidence of malignant tumors, it will be of great benefit when the number of such clinics in Japan and the rest of Asia is increased.

R3HDM1 haploinsufficiency is associated with mild intellectual disability.

R3HDM1 (R3H domain containing 1) is an uncharacterized RNA-binding protein that is highly expressed in the human cerebral cortex. We report the first case of a 12-year-old Japanese male with haploinsufficiency of R3HDM1. He presented with mild intellectual disability (ID) and developmental delay. He had a pericentric inversion of 46,XY,inv(2)(p16.1q21.3)dn with breakpoints in intron 19 of R3HDM1 (2q21.3) and the intergenic region (2p16.1). The R3HDM1 levels in his lymphoblastoid cells were reduced to approximately half that of the healthy controls. However, the expression of MIR128-1, in intron 18 of R3HDM1, was not affected via the pericentric inversion. Knockdown of R3HDM1 in mouse embryonic hippocampal neurons suppressed dendritic growth and branching. Notably, the Database of Genomic Variants reported the case of a healthy control with a 488-kb deletion that included both R3HDM1 and MIR128-1. miR-128 has been reported to inhibit dendritic growth and branching in mouse brain neurons, which directly opposes the novel functions of R3HDM1. These findings suggest that deleting both R3HDM1 and MIR128-1 alleviates the symptoms of the disease caused by loss-of-function mutations in R3HDM1 only. Thus, haploinsufficiency of R3HDM1 in the patient may be the cause of the mild ID due to the genetic imbalance between R3HDM1 and MIR128-1.

Clinical delineation, sex differences, and genotype-phenotype correlation in pathogenic KDM6A variants causing X-linked Kabuki syndrome type 2.

PURPOSE: The variant spectrum and the phenotype of X-linked Kabuki syndrome type 2 (KS2) are poorly understood. METHODS: Genetic and clinical details of new and published individuals with pathogenic KDM6A variants were compiled and analyzed. RESULTS: Sixty-one distinct pathogenic KDM6A variants (50 truncating, 11 missense) from 80 patients (34 males, 46 females) were identified. Missense variants clustered in the TRP 2, 3, 7 and Jmj-C domains. Truncating variants were significantly more likely to be de novo. Thirteen individuals had maternally inherited variants and one had a paternally inherited variant. Neonatal feeding difficulties, hypoglycemia, postnatal growth retardation, poor weight gain, motor delay, intellectual disability (ID), microcephaly, congenital heart anomalies, palate defects, renal malformations, strabismus, hearing loss, recurrent infections, hyperinsulinism, seizures, joint hypermobility, and gastroesophageal reflux were frequent clinical findings. Facial features of over a third of patients were not typical for KS. Males were significantly more likely to be born prematurely, have shorter stature, and severe developmental delay/ID. CONCLUSION: We expand the KDM6A variant spectrum and delineate the KS2 phenotype. We demonstrate that the variability of the KS2 phenotypic depends on sex and the variant type. We also highlight the overlaps and differences between the phenotypes of KS2 and KS1.

CDK19-related disorder results from both loss-of-function and gain-of-function de novo missense variants.

PURPOSE: To expand the recent description of a new neurodevelopmental syndrome related to alterations in CDK19. METHODS: Individuals were identified through international collaboration. Functional studies included autophosphorylation assays for CDK19 Gly28Arg and Tyr32His variants and in vivo zebrafish assays of the CDK19G28R and CDK19Y32H. RESULTS: We describe 11 unrelated individuals (age range: 9 months to 14 years) with de novo missense variants mapped to the kinase domain of CDK19, including two recurrent changes at residues Tyr32 and Gly28. In vitro autophosphorylation and substrate phosphorylation assays revealed that kinase activity of protein was lower for p.Gly28Arg and higher for p.Tyr32His substitutions compared with that of the wild-type protein. Injection of CDK19 messenger RNA (mRNA) with either the Tyr32His or the Gly28Arg variants using in vivo zebrafish model significantly increased fraction of embryos with morphological abnormalities. Overall, the phenotype of the now 14 individuals with CDK19-related disorder includes universal developmental delay and facial dysmorphism, hypotonia (79%), seizures (64%), ophthalmologic anomalies (64%), and autism/autistic traits (56%). CONCLUSION: CDK19 de novo missense variants are responsible for a novel neurodevelopmental disorder. Both kinase assay and zebrafish experiments showed that the pathogenetic mechanism may be more diverse than previously thought.

A novel missense variant in CUL3 shows altered binding ability to BTB-adaptor proteins leading to diverse phenotypes of CUL3-related disorders.

CUL3 forms Cullin-Ring ubiquitin ligases (CRL) with Ring-box protein and BTB-adaptor proteins. A variety of BTB-adaptor proteins have been reported to interact with the N-terminus of CUL3, which makes it possible to recognize various substrates for degradation. Regarding the association of CUL3 with neurodevelopmental disorders, a recent study reported three patients with global developmental delay, who carried de novo variants in CUL3. Here, we describe a novel de novo CUL3 variant (c.158G > A, p.Ser53Asn) identified in a patient with global developmental delay, who presented some novel dysmorphic features, including macrocephaly, characteristic facial features, and cutis marmorata. Immunoprecipitation and immunoblot analyses identified significantly weaker binding ability to some BTB proteins in CUL3-S53N compared to wild-type. Interestingly, label-free quantification proteomics analysis of samples immunoprecipitated by CUL3-S53N showed a significantly decreased interaction with some BTB proteins, while almost equal interaction or significantly increased interaction was observed with other BTB proteins. The binding between CUL3 and BTB proteins is essential for CRL substrate recognition, and alteration of their interaction is thought to result in the quantitative alteration in substrate proteins. It is possible that the difference of dysmorphic features between the present case and previously reported cases is caused by the distinctive effect of each CUL3 variant on substrate proteins. The clinical information of the present case will expand the picture of CUL3-related global developmental disorders, and subsequent cell biological analysis of the novel mutation will provide insight into the underlying molecular mechanism of how CUL3 pathogenic variants cause neurological disorders.

Role of chimeric transcript formation in the pathogenesis of birth defects.

Chimeric transcripts are formed by chromosomal aberrations. Little is known about the role of chimeric transcripts in the pathogenesis of birth defects. We reanalyzed RNA-seq data in alignment map files from the peripheral blood of 56 patients in whom the diagnoses could not be confirmed by standard exome analysis and transcriptome analysis to screen for chimeric transcripts using a dedicated software, ChimPipe. Chimeric analysis led to a diagnosis in two of the 56 patients: (a) the first patient had a chimeric transcript spanning the causative gene ZEB2 and the GTDC1 gene in its neighboring locus. RNA-seq revealed reads spanning exon 5 of ZEB2 and exon 7 of GTDC1. Whole genome sequencing revealed a 436-kb deletion spanning intron 4 of ZEB2 and intron 7 of GTDC1 and the diagnosis of Mowat-Wilson syndrome was made. (b) The second patient had a chimeric transcript spanning the causative gene KCNK9 and the TRAPPC9 gene in its neighboring locus. RNA-seq revealed reads spanning exon 21 of TRAPPC9 and exon 1 of KCNK9. Whole genome sequencing revealed a 186-kb deletion spanning intron 20 of TRAPPC9 and intron 1 of KCNK9 in this patient. KCNK9 gene is a maternally expressed imprinted gene. The diagnosis of Birk-Barel syndrome was made. Thus, both patients had chimeric transcripts that were directly involved in the pathogenesis of the birth defects. The approach reported herein, of detecting chimeric transcripts from RNA-seq data, is unique in that the approach does not rely on any prior information on the presence of genomic deletion.

Expanding the phenotype of biallelic loss-of-function variants in the NSUN2 gene: Description of four individuals with juvenile cataract, chronic nephritis, or brain anomaly as novel complications.

The NSUN2 gene encodes a tRNA cytosine methyltransferase that functions in the maturation of leucyl tRNA (Leu) (CAA) precursors, which is crucial for the anticodon-codon pairing and correct translation of mRNA. Biallelic loss of function variants in NSUN2 are known to cause moderate to severe intellectual disability. Microcephaly, postnatal growth retardation, and dysmorphic facial features are common complications in this genetic disorder, and delayed puberty is occasionally observed. Here, we report four individuals, two sets of siblings, with biallelic loss-of-function variants in the NSUN2 gene. The first set of siblings have compound heterozygous frameshift variants: c.546_547insCT, p.Met183Leufs*13; c.1583del, p.Pro528Hisfs*19, and the other siblings carry a homozygous frameshift variant: c.1269dup, p.Val424Cysfs*14. In addition to previously reported clinical features, the first set of siblings showed novel complications of juvenile cataract and chronic nephritis. The other siblings showed hypomyelination and simplified gyral pattern in neuroimaging. NSUN2-related intellectual disability is a very rare condition, and less than 20 cases have been reported previously. Juvenile cataract, chronic nephritis, and brain anomaly shown in the present patients have not been previously described. Our report suggests clinical diversity of NSUN2-related intellectual disability.

Biallelic loss of OTUD7A causes severe muscular hypotonia, intellectual disability, and seizures.

The heterozygous deletion of 15q13.3 is a recurrently observed microdeletion syndrome associated with a relatively mild phenotype including learning disability and language impairment. In contrast, the homozygous deletion of 15q13.3 is extremely rare and is associated with a much severer phenotype that includes epileptic encephalopathy, profound intellectual disability, and hypotonia. Which of the genes within the deleted interval is responsible for the more severe features when biallelically deleted is currently unknown. Here, we report a patient with profound hypotonia, severe intellectual disability, and seizures who had biallelic loss-of-function variants in OTUD7A: a 15q13.3 deletion including the OTUD7A locus, and a frameshift OTUD7A variant c.1125del, p.(Glu375Aspfs*11). Unexpectedly, both aberrations occurred de novo. Our experiment using Caenorhabditis elegans showed that worms carrying a corresponding homozygous variant in the homolog OTUB-2 exhibited weakened muscle contraction suggestive of aberrant neuromuscular transmission. We concluded that the biallelic complete loss of OTUD7A in humans represents a presumably new autosomal recessive disorder characterized by profound hypotonia, severe intellectual disability, and seizures.